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Abstract 

We study cutoff effects at tree-level of perturbation theory for maximally 
twisted mass Wilson, overlap and the recently proposed Creutz fermions. 
We demonstrate that all three kinds of lattice fermions exhibit the expected 
0(a 2 ) scaling behaviour in the lattice spacing. In addition, the sizes of 
these cutoff effects are comparable for the three kinds of lattice fermions 
considered here. Furthermore, we analyze situations when twisted mass 
fermions are not exactly at maximal twist and when overlap fermions are 
studied in comparison to twisted mass fermions when the quark masses are 
not matched. 



1 Introduction 



This paper is an investigation of the scaling behaviour towards the continuum 
limit for different kinds of lattice fermions at tree-level of perturbation theory. In 
a previous conference contribution [1] we had only discussed the case of Wilson 
twisted mass fermions [21 E]. See ref. [1] for a review on twisted mass lattice 
QCD. Here we add overlap fermions [5]. See ref. [6] for a review on chirally 
symmetric lattice actions and related topics. During the completion of this paper 
a new kind of lattice fermion appeared in the literature, the so-called Creutz 
fermions [7] and we decided to investigate also these lattice fermions, as well as 
the fermions defined by a related action suggested by Borici [8]. We will in the 
following generically refer to Creutz fermions having in mind both Creutz and 
Borici formulations. 

All of these three kinds of lattice fermions are expected to show an 0(a 2 ) scal- 
ing behaviour in the lattice spacing. While for overlap fermions this is achieved 
by an exact lattice chiral symmetry, for twisted mass fermions this is achieved by 
a residual Nf — 2 flavoured continuum chiral symmetry which needs, however, 
a tuning of the (untwisted) PCAC quark mass to zero. Creutz fermions exhibit 
an exact Nf — 2 flavoured continuum chiral symmetry, but they break however a 
number of discrete symmetries such as parity, charge conjugation and time reflec- 
tion [9]. This could lead to a plethora of terms in the Symanzik effective action, 
which would make approaching the continuum limit a rather difficult task. 

Wilson twisted mass fermions are described by an ultra-local action with only 
nearest neighbour interactions. Hence, they are comparably cheap to simulate. 
Their major drawback is the explicit 0(a 2 ) isospin-breaking. Overlap fermions 
have the great advantage of an exact lattice chiral symmetry [10J. However, 
while still local in the field theoretical sense [TT], they couple each lattice point 
with all others and are an order of magnitude more expensive to simulate than 
twisted mass fermions. Finally, Creutz fermions [7] are also described by an 
ultra-local action, again with only nearest neighbour interactions. As has been 
shown recently, they break however a number of discrete symmetries and isospin 
symmetry [9]. 

Twisted mass and Creutz fermions are defined for Nf = 2 flavours of quarks. 
These two quark flavours can either be taken as mass degenerate, or, some explicit 
flavour breaking term has to be added. In contrast, overlap fermions can be used 
for a single quark flavour. 

The focus in this paper is to study the scaling behaviour with the lattice 
spacing. In particular, we are interested in the relative size of the cutoff effects 
comparing maximally twisted mass, overlap and Creutz lattice fermions. 

In particular, we will consider the scaling behaviour of the meson correlation 
function the pseudoscalar mass and decay constant as well as the nucleon mass. 
It is important to remark here that, at tree-level of perturbation theory there is, 
of course, nothing like a real meson or baryon. However, since we use the same 
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interpolating fields as in full QCD, we will take the freedom to use the notation 
of mesons and baryons throughout this letter. 

We will also address the question of the size of cutoff effects when twisted mass 
fermions are not exactly tuned to maximal twist. In addition, we investigate 
ratios of mesonic quantities built from twisted and overlap fermions when the 
quark masses are unmatched. 

While actual practical calculations with overlap and twisted mass fermions 
are already rather advanced, the formulation of chiral invariant fermions following 
Creutz is still very new. 

In [TJ we gave a first account of the size of cut-off effects at tree-level of 
perturbation theory for twisted mass fermions. Here we extend the techniques 
described in [T], following [12], to the case of overlap and Creutz fermions. A 
more detailed account of our earlier calculations and a pedagogical introduction 
to the techniques we have used, can be found in [13]. 



2 Lattice propagators 

In this section we provide the momentum space propagators for the different 
kinds of lattice fermions we have considered. They are the building blocks for the 
computation of the meson and nucleon correlation functions from which in turn 
the masses and decay constants are extracted. To fix the notation our setup is a 
hypercubic lattice of size L 3 x T with spacing a. 



Wilson twisted mass fermions 



The expression for the Wilson twisted mass (Wtm) fermion propagator in the 
twisted basis, at tree-level of perturbation theory (PT) and in momentum space 
can be derived from the Wtm operator given for example in ref. [4], and it is 
given by 

o i \ -W?/* 1 / + M(p)U f - a/ig 75^3 m 
E M Vl + M{pf + fl 2 q 

where 

p M = - sin(ap M ), p^ = - sin(^), M(p) = m + °- ^ P% ( 2 ) 

and where 1 and 1/ are the identity matrices in Dirac and flavour space. The 
structure in colour space has not been written since it is just an identity matrix 
at tree level of PT. The parameters m and \i q represent the untwisted and 
twisted quark masses, respectively. Maximal twist -in the case of tree-level of 
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perturbation theory- is achieved by setting m = 0. We then expect to have only 
0(a 2 ) lattice spacing effects in physical correlation functions [3]. 

Overlap fermions 

The expression for the overlap propagator in momentum space at tree-level of 
perturbation theory can be derived from the expression of the overlap operator 
given for example in ref. pU], and it is 

ov[P) {-t-^Ym-'Y^pl + Mipy u 

where: 
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p<f 




('-t)w1i-tE<!) (6) 



and 1 is the identity matrix in Dirac space. Note that in the case of overlap 
fermions we only discuss one flavour. Due to the existence of an exact lattice 
chiral symmetry, we again expect an 0(a 2 ) scaling behaviour towards the contin- 
uum limit, if the correlation functions are computed with the proper improved 
operators (see for example ref. |14|). 

Creutz fermions 

The Creutz-Dirac operator can be found in different forms in ref. [7], [HE]. We have 
rescaled the Creutz-Dirac operator with a factor R that we leave unspecified for 
the moment. As will be discussed below and in the appendix, this normalization 
factor R is needed to obtain at tree-level the correct continuum limit for the 
correlation functions we have studied. In the appendix we define all the relevant 
functions and we show that the Creutz operator can be brought to the form 

D c{p) =i 2^PM^/« + m ol, (6) 

where and are linear combinations of gamma matrices defined in the ap- 
pendix. We recall that Creutz fermions, as explicitly shown in the appendix, 
depend in general on two paramenters R and C . From our final expression of the 
Creutz operator (jHJ), we can obtain the quark propagator for Creutz fermions. In 
the continuum limit, the quark propagator can be written as 

sr\P) - ^ — - - + , m 

E M E P PuPn a P^ a Pn + E P PpPv a Pn a P» + m o 
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which is not the continuum Dirac propagator unless the values of C and R, 
whose dependence is hidden in the matrices a and 7^, are chosen properly. 

From equation (j2J) it is also clear that in order to obtain the correct continuum 
limit of the quark propagator, we have to impose that 7 = a T 7 is again a gamma 
matrix, which is equivalent to say that the following relation holds true 

{(o T 7)a«, {a T l)v} = °" av 7a} = 2 ^ &pp a P v = 2 S^. (8) 

p,a p 

Taking into account the form of the matrix a defined in the appendix, 7 is a 
gamma matrix if and only if the following two conditions hold 

^2 a-pp, a pu = for fi ^ v =>■ — = ±1 C 

p 

^ a Pili a pu = 1 for fi = v =>• R 2 

p 

Therefore, from this discussion, it can be concluded that the right continuum 
limit of the quark propagator can be obtained only when C = 3/y/lO and R = 2. 

However, motivated by the first version of [7], we have decided to study 
as well the case C = 3/VT4. In this case we know that 7 is n °t a gamma 
matrix, since the condition given by equation ((HJ) does not hold. R can be still 
determined and must be R = 2\/2 in order to satisfy equation (fTUl) . As expected, 
the analytical expression obtained for the quark propagator in this case does not 
correspond to the continuum one. 

But, as it will be shown in the next sections where we present our results, 
the continuum limits for the pseudoscalar correlation function, mass and decay 
constant turn out to be the correct ones for the two values of C studied here. 

There are two reasons to explain this behaviour. On one hand, in the cor- 
relation function only the contribution from the pole survives in the sum over 
the momenta. This implies that all the possible crossed terms of the momen- 
tum which appear in the quark propagator, whose origin is the non-realization of 
equation (jHJ), are cancelled. Moreover the properly chosen value of R normalizes 
the quark propagator in order to obtain the 'would-be' correct continuum limit 
in absence of crossed terms. 

Borici fermions 

The action suggested by Borici [8] is a slight modification of the Creutz-Dirac 
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operator and reads 



f b(p) = J] i7A - ^ J] + m 1, (11) 
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where are linear combinations of gamma matrices defined in the appendix. 
The corresponding propagator can be written as 

Sb(p) = E,W + ™8 • (12) 

where the functions G^ap), defined in the appendix, are trigonometric functions 
of their argument. 



3 Correlation functions 

In this section we give the expressions for the pseudoscalar and proton corre- 
lation functions at tree-level of perturbation theory. While we evaluated the 
pseudoscalar correlation function for all lattice fermions considered in this work, 
we use the proton correlation function only to demonstrate how the 0(a) im- 
provement works when the Wilson average or equivalently the mass average is 
performed in case standard Wilson fermions are considered. 

The interpolating fields describing the charged pions, 7i + and n~ are 

V ± (x) =V 1 (x)TiV 2 (x) (13) 

where V a [x) = , ip(x) n f5 I !rip(x), with a = 1,2,3, is the pseudoscalar density and 
r a are the standard Pauli matrices. The computation for Wilson twisted mass is 
performed in the twisted basis, where the form of the local operator ([TBI stays 
unchanged [2]. 

The quark propagator can be decomposed in terms of the gamma matrices as 

s( P ) = s u (p)i + J2s,(ph, (14) 

ft 

in the case of overlap and Creutz fermions, while for twisted mass fermions an 
additional term proportional to 75 is present 

s(p) = Su(p)i + s Mi» + s 5 (pH. (15) 

With such a decomposition, the pseudoscalar correlation function can be written 

c ® ^EEEE eMf5 >>^>>i). ( 16 ) 

Pi p' 4 P A 4 

x The Wtm propagator has also a flavour structure. For the pseudoscalar correlation function 
only a single flavour component is needed. 
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with fi = U, 1, 2, 3, 4 or fi = U, 1,2,3,4, 5 depending on the kind of fermion 
that is being considered. N c is the number of colours and N d is the number of 
Dirac components. T = aN^ and N± is the number of lattice points in the time 
direction. The expression in eq. ([TBI) can be evaluated as it stands, or, a time- 
momentum representation of the quark propagator can be obtained for a lattice 
with infinite time extension, by performing the integration over analytically, 
see [13] for a discussion in case of twisted mass fermions. We have mostly used 
the representation of eq. (TT6]) , but checked for twisted mass fermions the results 
also in the time-momentum representation. In the appendix we give the explicit 
expression for the Wilson twisted mass quark propagator in the time-momentum 
representation for a lattice with infinite time extent. 

The local interpolating field describing the proton is given b)H 

&> a (x) = -V2e abc [F a {x] C-SVaO] u a , c (x). (17) 

The expression for the time dependence of the proton correlation function for 
standard Wilson fermions is then 

Cp#{t) = ^^{Lu&ftt) +L 4 (p,q,t)}, (18) 

v <f 

with the definitions 

L v (p,q,t) = Su{-{p + q),t){{N d + l) Su{p,t)Su{q,t) + 

4 

+ (N d + 3)J2S,(p,t)S,(q,t)} (19) 



L^{p,q,t) = S fl (-(p + q),t){(N d + 3)S u (p,t)S u (q,t) + 

4 

+ (Ni + V^S^QSpfat)}. (20) 
n=\ 

Su and Sfj, are the components of the quark propagator for standard Wilson 
fermions, which can be obtained from the expression of the quark propagator for 
Wilson twisted mass fermions, given in the appendix in eq. (1241 ). just by setting 
fig = 0. These results have been cross-checked with a standard inverter on a cold 
gauge configuration. 



2 The Greek (Latin) letters denote Dirac (colour) components and u, d denote the flavour 
content. C is the charge conjugation matrix and [] denotes spin trace. 
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4 Scaling tests on correlation functions, masses 
and decay constants 

In this section we give some results from a scaling analysis in the lattice spacing. 
At tree-level a dimensionless quantity can be only a function of mL, a/L and am, 
where m here indicates generically the quark mass. To perform the continuum 
limit one can fix mL to a certain value and the remaining dependence of the 
dimensionless quantity will be then in a/L. The continuum limit is then obtained 
sending N = L/a to infinity. In the following we will set a = 1 and the 1/N 
and 1/N 2 dependence of the dimensionless quantities under investigation will 
correspond to 0(a) and 0(a 2 ) scaling violations. We remind that if not infinity 
the time extent T will be always set to be proportianal to L. 

In particular, we consider the correlation function C at a fixed physical time 
t/N, the pseudoscalar decay constant fps and the pseudoscalar and proton masses 
M. This leads us to consider the dimensionless quantities N 3 C(t/N), NM and 
Nf PS . 

We will start our discussion with an explicit demonstration of 0(a) improve- 
ment for standard Wilson fermions from the Wilson and mass averaging proce- 
dures. We then turn over to the comparison of the size of scaling violations from 
all three lattice fermions considered here. In particular, we will compare the 
meson correlation function as well as the corresponding meson mass and meson 
decay constant. 

4.1 Wilson average and mass average for standard Wilson 
fermions 

In Ref. [3] it has been demonstrated that when averaging physical observables 
computed with standard Wilson actions, i.e. setting fi q = in eq. having 
opposite signs of the quark mass m (MA) or opposite signs of the Wilson pa- 
rameter r (WA), these quantities are 0(a) improved^. Since WA and MA at 
tree-level are equivalent, we will in the following only discuss the MA. As the 
physical observable we consider the proton mass as it can be obtained as the 
effective mass from the correlation function in eq. (fTHl) using timeslices at t = 4N 
and t = AN + 1. For the computation we fix |A^m | = 0.8. In the left graph of 
fig. Q] the behaviour of the proton mass NM P as a function of is given when 
two standard Wilson regularizations, differing only in the sign of the quark mass, 
are used. 

The behaviour of the proton mass in fig. [T] is linear in 1/N showing the 
expected 0(a) scaling violations of standard Wilson fermions. The proton mass 
reaches, in both cases, its continuum value of NM P = 3|7Vm |, as expected in 

3 The MA is actually done taking into account the chirality of the correlation function under 
investigation [3]. This is practically irrelevant for our study. 
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the absence of interaction. However, the slopes with which the continuum value 
is reached i.e. the coefficients of the 0(a) cutoff effects differ in sign. The right 
graph of fig. ED shows the continuum approach of the proton mass when the proton 
mass is averaged over positive and negative values of the quark mass Nm . In 
this case, the proton mass is plotted against 1/N 2 and the linearity in 1/N 2 nicely 
shows the 0(a)-improvement when the MA is applied. 

In order to be more quantitative, we have used the following fitting functions 
to describe the analytically computed values of the proton mass: 

Vi = a + a 1 - + a 2 — , y 2 = b + b 1 — + b 2 — . (21) 

Here y\ (y 2 ) is the physical observable under consideration and its value in the 
continuum limit is given by the coefficient a (bo). We use two functional forms, 
the first formula of equation (l2Ti ) for a leading behaviour (standard Wilson 
fermions) and the second formula for 0(a)-improved quantities. 

The two lines in the left graph of fig. Q] originate from a fit to equation ([2Til and 
correspond to the proton mass obtained from the same Wilson actions differing 
only in the sign of the quark mass. From the plot it is clear that in both cases 
the value of the proton mass in the continuum limit is the same and the expected 
one at tree-level of PT. 

From the fit, the corresponding coefficients a\ turn out to be the same in 
magnitude but have opposite signs for positive and negative quark masses. Thus, 
performing the (MA), it is to be expected that the 0(a) effects cancel and the 
scaling behaviour changes from a to a 4j behaviour. This can indeed be seen 
in the right graph of fig. [TJ Inspecting the fit coefficients a 2 and 6 ls we find 
a 2 « b\ ~ 0.5. Therefore, the magnitude of the leading order cutoff effects does 
not only change from an 0(a) to an 0(a 2 ) behaviour but also the 0(a 2 ) effects do 
not increase when performing the Wilson average with respect to the standard 
case. 

One interesting observation is the mass dependence of the coefficients a\ and 
bi of eqs. ([2T]) for the pion and proton masses. We observe that in the case of 
standard Wilson fermions, the coefficient ai/a , which determines the relative 
size of the 0(a) cutoff effects, vanishes in the chiral limit proportionally to Nuiq. 
For Wtm fermions at maximal twist the coefficient b\/bo, which determines the 
relative size of the 0(a 2 ) cutoff effects, vanishes in the chiral limit, proportionally 
to (iV/i q ) 2 . 



4.2 Comparing maximally twisted mass, overlap and 
Creutz fermions 

One interesting question is the relative size of cutoff effects when comparing max- 
imally twisted mass, overlap and Creutz fermions. We have therefore performed 
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1/N 1/N 2 



Figure 1: In the left graph, the cutoff effects and the continuum limit of the proton mass 
obtained from two standard Wilson actions differing only in the sign of the quark mass, \Nmo\ = 
0.8 are shown. The lattices are 4 < N < 20. The lines are fits according to eq. IpTL In the 
right graph, the average of the proton masses obtained from the same two standard Wilson 
regularizations with quark masses Nttiq = ±0.8 (MA) has been calculated. The solid line 
represents a fit to eq. l(2T|). 



NM ps 


a 


b 


c 


MTM 


1 


-0.0104167 


0.000292154 


OVERLAP 


1 


0.0208333 


0.000783943 


BORICI 


1 


-0.0494786 


0.00558893 


CREUTZ - VlO 


1 


-0.00781168 


-0.010171 


CREUTZ - y/U 


1 


-0.0488288 


0.00287405 



Table 1: Table of fit coefficients for the pseudo scalar mass using NMps — a + bj^ + Cjjz- 



a scaling analysis for the correlation functions themselves at a fixed physical dis- 
tance, the pseudoscalar mass and the pseudoscalar decay constant. Since all the 
quantities under investigation are 0(a)-improved, we show them all as a function 
of 1/N 2 . 

Let us start with the pseudoscalar mass in fig. [3 As expected, we indeed find 
a nice linear behaviour of the mass as a function of 1/N 2 . However, we observe 
that the different kind of lattice fermions we have studied show quite distinct 
lattice artefacts at 0(a 2 ). This can also be seen in fig. where we exhibit the 
dependence of the correlation function and the pseudoscalar decay constant as a 
function of 1/N 2 . 

In order to see the cutoff effects in a quantitative way, the coefficients of the fit 
and the corresponding fit function for the three lattice quantities here considered, 
obtained from the different regularizations are presented in Table [U Table [2] and 
Table El respectively. 

From the results obtained here, no clear picture of a particularly good or bad 
fermion discretization emerges. While we find that indeed all three kind of lattice 
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Figure 2: The cutoff effects and the continuum limit of the pseudo scalar mass is shown. We 
compare overlap, maximally twisted mass and Creutz fermions at fixed quark mass, Nm = 0.5, 
N i^L q = 0.5 and NrriQ = 0.5, respectively. The lattices are 4 < N < 64. 
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Figure 3: In the left graph, the cutoff effects and the continuum limit of the pseudo scalar cor- 
relation function is shown. We compare overlap, maximally twisted mass and Creutz fermions 
at fixed quark mass, Nm = 0.5, N^i q = 0.5 and Nmo = 0.5, respectively. The lattices are 
4 < N < 64. In the right graph, the pseudo scalar decay constants obtained from the same 
lattice fermions as in the left graph are shown. 



Nfps 


a 


b 


c 


MTM 


3.4641 


0.0541266 


-0.000815548 


OVERLAP 


3.4641 


0.108253 


0.00554908 


BORICI 


3.4641 


-0.0676637 


-0.00486739 


CREUTZ - 


3.4641 


0.293217 


-0.0770494 


CREUTZ - VlA 


3.4641 


-0.00790885 


-0.0367598 



Table 2: Table of fit coefficients for the pseudo scalar decay constant using jV/ps = a + bjj? + 
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N' A C(t/N = 4) 


a 


b 


c 


MTM 


0.109894 


0.00457891 


-3.33302 ■ 10~ 5 


OVERLAP 


0.109894 


0.0045789 


0.000181822 


BORICI 


0.109894 


0.00114427 


-0.00135812 


CREUTZ - VlO 


0.109894 


0.0194625 


-0.00286602 


CREUTZ - y/U 


0.109894 


0.00486428 


-0.002942 



Table 3: Table of fit coefficient for the pseudo scalar correlation function using N 3 C = 



fermions show the expected 0(a)-improvement, the (relative) size of the 0(a 2 ) 
effects depends pretty much on the observable considered. If at all, one could say 
that maximally twisted mass fermions show uniformly small 0(a 2 ) cutoff effects. 
On the other hand, it is somewhat amazing that Creutz fermions which break 
a number of important discrete symmetries do not suffer from very large 0(a 2 ) 
cutoff effects. From our scaling analysis it is not possible to exclude a certain 
type of lattice fermion. Only scaling tests for the interacting theory will reveal 
the size of actual scaling violations of the observable considered. 

5 Effects from non-optimal tuning 

This section is devoted to the question of effects when tuning is performed non- 
optimally. In particular, we study the cutoff effects when there is an O(a) error 
in tuning to maximal twist. As a second example, we consider the case when the 
quark masses of two lattice fermion formulations are not exactly matched. This 
case is relevant for so-called mixed action simulations. 

5.1 Out of maximal twist 

Here we want to study a situation when we allow an O(a) error in setting the 
untwisted quark mass to zero. In order to realize this situation at tree-level of 
perturbation theory we 'force' these effects by simply fixing the twisted mass to 
be the physical quark mass and the untwisted mass is set to be proportional to 
as N/j, q = a and Nm — 4 ^ 0(a) where a is kept fixed and (3 is a measure 
parametrizing the amount of violation of the maximal twist setup. The twist 
angle u> and the bare polar mass M can be obtained as a function of a and (5 as 

Therefore, even if the condition of maximal twist can only be obtained up to O(a) 
cutoff effects, which is generally the case in practical numerical simulations, the 
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0.03 
1/N 2 



0.001 
1/N 2 



Figure 4: Left graph: Behaviour of the pion mass as a function of -j^, for lattices with size 
4 < N < 64. The twisted quark mass is set to Nfi q = 1.0 and the untwisted quark mass is zero 
up to O(a) cutoff effects i.e. Nm = jj with /3 = 0.0, 1.0, 2.0, 10.0. Right graph: a zoom of the 
graph on the left with an additional fit for the analytical data corresponding to /? = 10.0 which 
considers only large lattices 40 < N < 64. 



observables, which are only functions of the polar mass, are still automatically 
O(a) improved. 

Moreover, equation ( 1221 also shows how the size of the leading discretization 
effects depends on the ratio between the untwisted and twisted quark masses. 
This ratio in turn determines the value of the lattice spacing at which the 
asymptotic scaling sets in. Only when this ratio is small enough and hence 
the lattice does not need to be chosen too large a reliable continuum limit using 
reasonably sized lattices can be performed. The left graph of fig. [4] demonstrates 
that the asymptotic scaling sets in for lattices with 4 < N < 64 when < £ < 2. 

However, for ^ > 10 the continuum limit is not reliable anymore if N is chosen 
to be too small. This can be seen in the right graph of fig. [H Using only small 
values of iV leads to an inconsistent continuum limit value. Therefore, larger 
lattices are needed in order to obtain the correct continuum behaviour as can be 
also seen in the right graph of fig. BJ Here we have added a fit of the data for a 
value of (3 = 10.0 taking only large lattices into account, i.e. using only values of 
N > 40. In this case, indeed the right continuum value is obtained. Of course, 
for practical simulations, using only lattices with N ^> 40 appears to be rather 
unrealistic so to chose a correct ratio - becomes important. 

5.2 Unmatched quark masses 

In this section we want to study the continuum limit and the size of the cutoff 
effects of lattice quantities constructed from ratios of physical observables com- 
puted on the lattice from two different regularizations i.e. here Wilson twisted 
mass fermions at maximal twist and overlap fermions. In particular, we want to 
study the situation when both quark masses are not exactly fixed to the same 
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0.8 1 1 1 1 ' 1 1 1 0.9995 1 1 1 1 1 1 

0.01 0.02 0.03 0.04 0.05 0.06 0.0001 0.0002 0.0003 0.0004 0.0005 

1/N 2 1/N 2 

Figure 5: The cutoff effects and continuum limit of the ratio of the pseudoscalar mass 
computed for maximally twisted mass and overlap fermions. In both graphs Nfi q = 0.5, 
Nm = 0.5 — 0.4/iV 2 and t/N = 4. The left graph shows the full range of lattice sizes considered 
while the right graph represents a zoom. 



value but differ up to 0(a 2 ) effects. 

The reason for studying such setup is that in real simulations using a mixed 
but 0(a)-improved action both masses can be fixed to the same value only up to 
0(a 2 )effects. In order to realize non-matched quark masses, we fix the twisted 
quark mass exactly at Nfi q = 0.5 and allow for an 0(a 2 ) error in setting the 
overlap quark mass, 

Nm = 0.5 -v/N 2 . (23) 
We will vary the parameter v from v — to v — 4.0. 

Similar to the case of twisted mass at a non-optimally tuned value of the 
bare quark mass, the leading 0(a 2 ) cutoff effects can become very large when the 
quark masses are substantially mis-matched. This might even induce the danger 
of not achieving the right continuum limit when only too small lattice sizes are 
used. 

This is exemplified in fig. [5l In the left graph we show the scaling behaviour 
in l/iV 2 for various values of the parameter v for the full range of the lattice 
sizes. Ifv is small, we observe 1/iV 2 scaling violations. Only in the case v = 4, do 
the 1/N 4 corrections become relevant. More quantitatively a linear fit in 1/N 2 
including all the data points starting fron N = 4 is reliable for v < 1. While for 
v ~ 4 we can include, in a linear fit in 1/N 2 , only the data points starting from 
N= 12. 

Fig. [5] demonstrates that not only the slope of the ratio of pseudoscalar masses 
versus 1/N 2 increases as we make the value of v larger, but also that the 0(a 4 ) 
effects become significant as can be seen from the curvature in the plot for v = 4. 
Moreover, for v = 4, the continuum limit fit is not reliable anymore if too small 
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lattices are considered, which can be seen in the right graph of fig. El Thus, if 
the matching is not performed with a good enough precision, large lattices have 
to be used to avoid higher order cutoff effects and perform a reliable continuum 
limit. 

6 Conclusions 

In this paper we have performed a scaling test in the lattice spacing towards the 
continuum limit for three kinds of lattice fermions. Our setup has been tree-level 
of pertubation theory. The lattice fermions considered here were Wilson twisted 
mass, overlap and Creutz and Borici fermions. We looked at the pseudoscalar 
correlation function at a fixed time and the corresponding pseudoscalar mass and 
decay constant. 

As a first step, we have verified automatic 0(a) improvement for Wilson 
twisted mass fermions and showed, with the example of the proton mass, the 
mechanisms of mass average. 

The relative comparison of all three kinds of lattice fermions we have consid- 
ered did not result in a clear picture in the sense that one lattice fermion shows 
consistently smaller or bigger 0(a 2 ) lattice artefacts than the other. Rather we 
found that the sizes of 0(a 2 ) lattice artefacts depend on the observable consid- 
ered with perhaps the exception of maximally twisted mass fermions which shows 
a rather uniform behaviour with small 0(a 2 ) effects. Therefore, we expect that 
also in practical simulations the 0(a 2 ) lattice artefacts can turn out to be quite 
different depending on the physical observable considered. 

Finally, we studied the situation when parameters are tuned non-optimally. 
We considered twisted mass fermions when an 0(a) error was allowed in the 
tuning to maximal twist. In addition, we looked at ratios of physical quantities 
built from overlap and twisted mass fermions. In this case we allowed an 0(a 2 ) 
error in the matching of the quark masses. Our conclusion of these studies is that 
when the corresponding 0(a 2 ) error is too large, the continuum limit becomes 
not reliable if the lattice spacing is not small enough. 
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Appendix 



In this appendix we summarize some of the formulae and definitions which have 
been used in this work. 

Wilson twisted mass fermions 

For Wilson twisted mass fermions we give the analytic expression for the time- 
momentum representation of the quark propagator on a lattice with infinite time 
extent. Details of the computations can be found in ref. [13j. 

The propagator reads^, 

S(p,t) = — — \sgn(t) sinh Eu^lf - ia)C(p)lf + [(1 - cosh Ex) + aM(p)] 11/ 

- ia/Vr5T 3} e - a iH +%0 _1_ 1 1 / 



(24) 



wit 



itlJl 



3 

a 



M(p) = m + - V\ 1C(p) = Yl 7iPi, ( 25 ) 



2 

i=l i=l 



C ° Sh El = 1 + 2(l + oM(p)) ' (26) 

£ 1 = \l + aM(p)\smhE 1 . (27) 



Creutz fermions 

Our starting point for the analysis of the Creutz action is the operator given 

by 

Dc(p) = i Y + c»( a P)) 7 M + m o 1 (28) 

where 1 is the identity matrix in Dirac space and the momentum has already 
been reexpressed in terms of the pole as 

?m = 9+2V- ( 29 ) 



4 This expresion is valid for all the possible values of the discrete Euclidean time t, even if 
it is negative or zero. This generalization has been done by using the properties of the Fourier 
transform under change of sign of the argument. 

5 sgn(i) is the sign of t, and we have denoted sgn(O) = 0. It is just a convention in order to 
give one general expression for the propagator for all possible values of t. 
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Starting from a reciprocal lattice labelled here by the poles of the Creutz-Dirac 
operator are localized at (q, q, q, q) and (—q, —q, —q, —q), where q is related to the 
parameter C by the equation C = cos(g). In the following, we consider only the 
former pole. In addition, we also define the parameter S = sin(g). 

The relevant trigonometric functions for Creutz fermions are defined by 



3S 



s 4 (ap) 



c k (ap) 



S 

RC 



c k (ap) 



RC 
3 

Ci(ap) = —c A {ap) 



with the functions s and c given by 



si(ap) 
s 2 {ap) 
s 3 (ap) 
s 4 (ap) 

ci(ap) 
c 2 (ap) 
c 3 (ap) 
c A (ap) 



[Pi +P2-P3- fa] 
[Pi -P2-P3+ Pi] 
[Pi -P2+P3- fa] 

[-fa -fa- fa- fa] 

2,-2 -2 -21 
~>2 ~ Ps - Pa] 

a 



a r "2 
-2^1 



1*2 "2 -2 1 "21 

[Pi -P2-P3+ Pa] 



" \ -2 '2,-2 -21 

~2 LPi -P2 +Ps -P4J 

a r-2 1 -2 1 -2 1 -21 
~2 LPi + ^2 + P3 + Pa\ 



(30) 
(31) 

(32) 
(33) 
(34) 
(35) 

(36) 
(37) 
(38) 
(39) 



With the expression of eq. ( l28l) we obtain the corresponding propagator with 
standard manipulations as 



Sc(p) 



-i E M [s^ap) + c M (ap)) 7^ + m 1 



(40) 



In order to understand the continuum limit of the Creutz-Dirac operator and the 
role of the factor R, we will use in the following the notation of Borici in [8], with 
a scalar product (7,2;) = 5^„7^%- 
We thus write eq. ( 1281 ) as 

i (a T 7, s(ap) + a" 1 bc(ap)) + m t (41) 



-Dc(p) = i (7, as(ap) + bc(ap)) + m 1 
where 



s(ap) = - 
a 



h,P2,P3,Pi) 



c(ap) = (Pi ^2^3:^4) 



2\ T 



(42) 



and the matrices a and 6 



1 

7? 



1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 


1 


-1 


3S 


3S 


35 


3S 


C* 


C 


C 


C 
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/ 1 



1_S 
RC 



1 \ 



V 



1-1-1 i 
l-ii-i 

3C 3C 3C 3C 

s s s s 



J 



As a next step, we define 

a = a -1 b , 7 = a T 7 , f = «7 = aa T 7 
from which we obtain the final form of the Creutz-Dirac operator 



(43) 



(44) 



This is the expression given in eq. ([6]). 



Borici fermions 

Here we give the definitions we have used for Borici fermions. The T matrices 
multiplying the Wilson term are defined by 



(45) 



with the definition 



1 

a= 2 



( 1 -1 -1 -1 \ 
-1 1 -1 -1 
-1 -1 1 -1 

V-i -l -l i j 

We incidentally note that X] m 7m = ~ S/x-'-V 

The trigonometric functions appearing in the Borici propagator are given by 



(46) 



Gi(ap) 


= Pi- 


a |~"2 i -2 -2 -21 

1 [Pi +P2-P3- Pa\ 


G 2 {ap) 


= P2~ 


a r -2 1 -2 -2 '2i 
4 [-Pi +P2-P3- Pa] 


G 3 (ap) 


= P3~ 


a r '2 "2,-2 "21 

4 [-Pi - P2 + Ps - P4J 


G A (ap) 


= p4 - 


fl r "2 "2 "2 1 "21 

1 [~Pi -P2-P3+ Pa\ 



(47) 
(48) 
(49) 
(50) 
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